Exclusive vector meson production at HERA from QCD with saturation 
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Following recent predictions that the geometric scaling properties of deep inelastic scattering 
data in inclusive "/*p collisions are expected also in exclusive diffractive processes, we investigate the 
diffractive production of vector mesons. Using analytic results in the framework of the BK equation 
at non-zero momentum transfer, we extend to the non-forward amplitude a QCD-inspired forward 
saturation model including charm, following the theoretical predictions for the momentum transfer 
dependence of the saturation scale. We obtain a good fit to the available HERA data and make 
| predictions for deeply virtual Compton scattering measurements. 

£ 

I. INTRODUCTION 

Geometric scaling Q is a striking feature of deep inelastic scattering (DIS) data at small Bjorken x, or at large 
, rapidity Y — log(l/x). The photon-proton total cross section (T^ ot p_> (x, Q 2 ) obeys a scaling in the single variable 
Q 2 /Q 2 (x) where Q is the virtuality of the photon and Q s the saturation scale. This momentum scale increases with 
t— ( \ rapidity according to Q 2 (x) — QqX~ x with A ~ 0.3 and Qo ~ 0.1 GeV (giving a saturation scale of order 1 GeV 
for x ~ 10 4 ). Recently, the diffractive cross-section aJ i ^~* Xp and the elastic vector-meson production cross-section 

Cn| ' a VM* VP were shown to also feature scaling behaviours, with the same saturation scale Q. 

From the theoretical point of view, the high-energy or small— x limit can be studied within the QCD dipole 
■ picture 0]. Introducing the dipole-proton elastic scattering amplitude T(r 7 Y), where r is the transverse size of 
the dipole, the corresponding law T(r,Y) = T(rQ s (Yj) appears to be a genuine consequence of saturation effects 
[3] characteristic of the high-density regime of perturbative QCD. More precisely, the evolution towards saturation is 
O conveniently described by the non-linear Balitsky-Kovchegov (BK) equation Q that resums QCD fan diagrams in the 
leading- logarithmic approximation. Geometric scaling for T(r,Y) appears to be a mathematical consequence of the 
asymptotic solution in Y of this nonlinear equation in terms of travelling waves Q . Within this theoretical framework, 
Q S (Y) — Qq£I s (Y) and the rapidity dependence of the saturation scale, given by il s (Y), is then obtained from features 
of the Balitsky-Fadin-Kuraev-Lipatov [7[ (BFKL) kernel, which drives the linear regime of the BK evolution. 

The total cross section u2 ^~* X (as well as the diffractive cross section o-J i ^~* Xp and the elastic vector-meson produc- 
tion cross section o^y[* Vp ) is related to the forward elastic dipole-proton amplitude. It is natural to ask the question 
whether geometric scaling survives when considering non- forward amplitudes. To that purpose, it is temptative to 
study the BK equation which is written for T(r,h;x), the dipole-proton elastic scattering amplitude, depending not 
only on the dipole size r but also on the impact parameter b. However, the aforementioned scaling law cannot be 
easily adapted to account for the impact parameter dependence. The analysis of the BK equation in impact-parameter 
space @, Q even shows a contradiction with confinement: the large-b dependence of the solutions develops a power-law 
tail with decreasing x. 

Previous studies [13, EH have shown that the travelling-wave method can be extended to the BK equation with 
transverse-space kinematical dependence provided one investigates the problem in terms of momentum transfer instead 
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FIG. 1: Vector-particle production in the dipole frame. The amplitude factorises as a product of vertex functions and elastic 
dipole-proton interaction amplitude (see formula (O). 



of impact parameter. The remarkable point [lfj is that the travelling-wave analysis of the BK equation can be better 
achieved in momentum space, by Fourier transforming T(r, b; Y) into T(r, q; Y) where q is the momentum transfer. 
From the knowledge of the exact solutions of the BFKL equation at non-zero momentum transfer , the travelling- 
wave property has been extended [l(J (for large Q): 

• at small momentum transfer, i.e. when |r||q| < |r|Qo < \r\Q ~ 1, one has asymptotically T(r,q;Y) = 
T(\r\Q n s (Y) , q) , recovering the forward result; 

• at intermediate, semi-hard, momentum transfer, i.e. when |r|Qo < |r||q| < |r|Q~ 1, one has asymptotically 
f(r,q; Y)=f (|r||q|n s (y),q); 

• at large, hard, momentum transfer, i.e. when |r|Qo < l r IQ < j cj| < — ' 1, one has no more saturation, (giving an 
explanation of the impact-parameter puzzle). 

This introduces a q-dependent saturation momentum in the near-forward and intermediate transfer region with a 
rapidity evolution il s (Y) keeping the same form as in the forward case. Those predictions were confirmed both by 
analytical and numerical analysis of the BK equation completely formulated in momentum space [llj |. 

Our purpose here is to investigate an interesting phenomenological prediction from these results: the geometric 
scaling property should manifest itself in exclusive vector meson production and deeply virtual Compton scattering 
(DVCS), which are experimentally measured at moderate non-zero momentum transfer t = — q 2 . In this paper, we 
analyse whether or not the available data from HERA are sensitive to a t— dependent saturation scale. To do this, we 
use a QCD-inspired saturation model for the dipole amplitude T(r, q; Y). Since it is important to include the charm 
in the analysis, both for its impact on DVCS, non-charmed mesons and also J production, we make use of a recent 
saturation model that successfully includes charm Note that our parametrisation, which uses the momentum 

transfer q instead of the impact parameter b as suggested by this derivation of saturation effects in perturbative 
QCD, provides a fruitful phenomenological framework since the data are directly measured as a function of t — — q 2 . 

The plan of the paper is as follows. In section II, we recall the formulation of vector meson production and DVCS 
differential cross-sections in terms of the dipole scattering amplitude T(r,q;x). In section III, we briefly explain the 
asymptotic travelling-wave properties of this amplitude and how they translate into geometric scaling at non zero 
momentum transfer; we also introduce our QCD-inspired model for T. In Section IV, we present fits to the vector- 
meson production experimental data, discuss our results and present predictions for DVCS. Section V concludes. 

II. EXCLUSIVE VECTOR MESON PRODUCTION IN DIS AT SMALL x 

In the small— a; limit, it is convenient to describe the scattering of the photon in a particular frame, called the dipole 
frame, in which the virtual photon undergoes the hadronic interaction via a fluctuation into a colorless qq pair, called 
dipole, which then interacts with the target proton. The wavefunctions ipj ^ ^(z, r; Q 2 ) describing the splitting of a 

virtual photon with polarization A = 0, ±1 into a dipole are well known. The indices h = ±l and h~±l denote the 
helicities of the quark and the antiquark composing the dipole of flavor /. The wavefunctions depend on the virtuality 
Q 2 , the fraction z of longitudinal momentum (with respect to the j*—p collision axis) carried by the quark, and the 
two-dimensional vector r whose modulus is the transverse size of the dipole. Explicit formulas for the QED functions 
^Jhh can ^ e f° un d i n the literature (l4j | and are recalled in Appendix lAl 



The exclusive production of vector mesons is represented in FigfT] the photon splits into a dipole of size r which 
scatters elastically off the proton, with momentum transfer q, and recombines into a vector meson whose mass we 
shall denote My. To describe this process, we need to introduce the wavefunctions i/iJ^(z,r;My) which describe 
the splitting of the vector meson with polarization A into the dipole. In fact, to compute the vector-meson production 



amplitude pictured in Fig[TJ we need the transverse (T) and longitudinal (L) overlap functions $ 
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These functions depend on the meson wavefunctions i/j , and different models exist in the literature [TH, [H, [l7| • We 
shall discuss two different choices later in this paper: the boosted Gaussian (BG) wavefunctions [l5j and the light-cone 
Gaussian (LCG) wavefunctions [r^| . For completeness, we give explicit expressions for those overlap functions in 
Appendix El 

If q denotes the transverse momentum transfered by the proton during the collision, the differential cross-section 
with respect to t = — q 2 reads 
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In that expression, A refers to the imaginary part of the scattering amplitude and is given by 
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(3) 



where x and y are respectively the transverse positions of the quark and the antiquark forming the dipole. T(x, y; Y) 
is the dipole-proton scattering amplitude and carries all the energy dependence via the rapidity Y which is now 
obtained from the centre-of-mass energy W and the vector-meson mass My using 



Y = log 



W 2 + Q< 
Ml + Q< 



The prefactor 1 + 1 in ([2]) accounts for the contribution coming from the real part of the amplitude and can be 
obtained using dispersion relations: 
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It is convenient to consider the dipole-proton scattering amplitude as a function of r = x— y and b = zx+(l — z)y 
and to introduce the following Fourier transform: 



f(r, q; Y) = J d 2 b e lqb T(r, b; Y) 



(5) 



Indeed, as we shall discuss in the next section, this quantity features the geometric scaling property at non-zero 
momentum transfer. The formula we shall use finally reads 
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dz $jJ(z,r;Q 2 ,M 2 ) e"^ r f(r,q;F) 



(6) 



Note that formula ((6]) can also be used to compute the DVCS cross-section a 1 provided one uses the vertex 

for a real photon instead of a vector meson in (|TJ) . The overlap function between the incoming virtual photon and the 
outgoing transversely-polarized real photon is now model-independent and given by 



r; Q 2 ) = W,hl^ r; °)] ^ ^/,a^> r; 
fhh 



(7) 



III. GEOMETRIC SCALING AT NON-ZERO MOMENTUM TRANSFER 



We have expressed the exclusive vector-meson production cross-sections © in the high-energy limit in terms of 
the Fourier-transformed dipole scattering amplitude off the proton T(r,q;Y). Its evolution towards large values of 
Y is computable from perturbative QCD and the most important result about the growth of the dipole amplitude 
towards the saturation regime is probably the geometric scaling regime. It first appeared in the context of the proton 
structure function, which involves the dipole scattering amplitude at zero momentum transfer. At small values of x, 
instead of being a function of a priori the two variables r = |r| and Y, the dipole scattering amplitude is actually a 
function of the single variable r 2 Q 2 (Y) up to inverse dipole sizes significantly larger than the saturation scale Q S (Y). 
More precisely, one can write 

f(r,q = 0;y) = 27ri?2jV( r 2 Q 2 (r)) > (8) 

implying (for massless quarks) the geometric scaling of the total cross-section at small x : a^ Q ^ X (Y,Q 2 ) — 

a tot^ X ( T = Q 2 IQ 2 s (y))- This has been confirmed by experimental data [l[ (see also for the geometric scaling 
of the diffractive and elastic vector-meson production cross-sections) with Q 2 (Y) — Q^e ', and the parameters 
A-0.3 and Q ~0.1 GeV. 

As explained in the Introduction, it has been shown in a recent work [ToL IllT ] that the geometric scaling property 
can be extended to the case of non zero momentum transfer, provided r|q| <C 1. We obtained that equation ((HI) can 
be generalised to 

f (r, q; Y) = 2nR 2 p /(q) N(r 2 Q 2 s (Y, q)) , (9) 

with the asymptotic behaviours Q 2 (Y, q) ~ max(Qo,q 2 ) exp(AV) and an unknown form factor /(q) of non- 
perturbative origin. 

In practice, we need to specify three ingredients: 

• The q dependence of the saturation momentum is parametrised as 

Q 2 s (Y,q) = Q 2 (l + cq 2 )e XY (10) 
in order to interpolate smoothly between the small and intermediate transfer regions. 

• The form factor /(q) catches the transfer dependence of the proton vertex. It has to be noticed that this form 
factor is factorised from the projectile vertices and thus does not spoil the geometric scaling properties. For 
simplicity, we use /(q) = exp(— Bq 2 ). 

• The scaling function N is obtained from the forward saturation model [HI, [l|| : 



iV i^Hy exp 



for r 2 Q 2 s (Y) < 4 



2kXY 

N(rQ s (Y),Y)= { ^ , (n) 

1 - e -^{^Ql(Y)) for r 2Q2 (y) > 4 

with a and (3 uniquely determined from the conditions that N and its derivative are continuous at r 2 Q 2 (x) = 4. 
The amplitude at the matching point is chosen to be No = 0.7. 

In this work, we shall consider the same form for the forward amplitude as in the Iancu-Itakura-Munier (IIM) 
saturation model [l8[ but in its version extended in [l3| to include heavy quarks (with m c = 1.4 GeV, = 
4.5 GeV, and m/ = 0.14 GeV for the light flavors). The coefficient k = 9.9 is obtained from the BFKL kernel 
while the critical exponent 7 C = 0.7376 is fitted to the HERA measurements of the proton structure function, 
along with the remaining parameters. The saturation scale parameters are A = 0.2197 and Qo = 0.298 GeV and 
the proton radius is i? p = 3.34 GeV" 1 . Note that, after including heavy quark contributions, the saturation scale 
stays above 1 GeV for x=10 -5 , rather than dropping to about 500 MeV as is the case in other studies (see the 
discussion in [131]). This parametrization is also successful in describing inclusive diffraction data [191 ]. 

The last factor in the expression for small dipole sizes introduces geometric scaling violations, important when Y 
is not large enough, as predicted by the high-energy QCD evolution. It controls the way how geometric scaling is 
approached. Indeed, it can be neglected for log(r 2 Q^/4) < y/2nXY, meaning that the geometric scaling window 
extends like VY above the saturation scale (in logarithmic units). 



The final expression of our QCD-based saturation model is thus 

f (r, q; x) = 2nR 2 p e~ B * ''N(rQ s (x, q), Y) (12) 

which is an extension of the forward model (T3 , fl8| including the QCD predictions for non zero momentum transfer. 
Indeed, this model reproduces the initial model for q = and ensures that the saturation scale has the correct 
asymptotic behaviours. We have two parameters: c related to the scale at which the q— dependence of the saturation 
scale becomes important and B, the t— slope of the form factor that we have taken as simple as possible. Those 
parameters have to be fitted to the experimental measurements of elastic vector-meson production as we shall comment 
in details in the next section. 

Note that the form factor depends on a single slope B, independently of the vector-meson in the final state. This is 
in agreement with predictions from the BK equation flO . [ll[ which implies the factorisation of the non-perturbative 
contribution in (|12p. The difference in the processes of p, <j> and J/^ production (or even DVCS cross-sections) fully 
comes from the vector-meson wavefunctions. Through their wavefunction, different vector mesons are sensitive to 
different dipole sizes, and thus feel saturation differently. 

IV. FIT RESULTS AND DISCUSSION 

In this section, we compare our formulation ([6]) and (|12j) for vector-meson production with the experimental 
measurements of p, (j> and J/\f r elastic production. We shall first state explicitly the data used to perform the fit, then 
give more details concerning the fit itself, present our results and compare them with other possible approaches, and 
finally we shall give predictions for DVCS measurements. 

A. Data selection 

In this analysis, we shall include the differential cross-section da/dt as well as the elastic cross-section <j c \. The 
dataset is obtained from the p-meson production measured by HI [20| . the 0-meson production measured by ZEUS 
[2lj . and the J/^-meson production measured by ZEUS [13] and HI |23l|. We have not taken into account the ZEUS 
data for p mesons [24| since they lie in the low-Q 2 region where we do not expect the condition r|g| ~ \J\AIQ 2 ^ 1 
to be valid. The ratio R — vl/vt between the longitudinal and transverse elastic cross-sections is mostly dependent 
of the choice of wavefunction, hence we have not included it in the analysis. Finally, this gives a total of 269 data 
points which have to be reproduced with the two parameters B and c in equation (|12j) . Indeed, we keep the same 
parameters as in fl3j for t — as they provide a good description of the inclusive F% data. 

B. Meson wavefunctions and alternative models 

In order to test the quality of our parametrisation with t— dependent saturation scale, we shall compare our results 
with two alternative parametrisations. Since the main aim of this work is to analyse the t— dependence of the 
saturation scale, the first alternative we shall consider is a model in which we only adjust the slope B while keeping 
the saturation scale independent of the momentum transfer. 

For the second alternative, we notice that if one assumes a slope B for the differential cross-sections at a given W 
and Q 2 , the slope is experimentally observed to depend on Q 2 . Within our parametrisation, it is the t dependence 
of the saturation scale which is expected to account for this behaviour in Q 2 . This will be compared with a model 
where we keep Q s constant in t (c = 0) and explicitly introduce a Q 2 dependence of the slope B. We shall use 
B(Q 2 ) = B + B' /(Q 2 + My), which shows the same trend as the experimental estimations of B(Q 2 ). Note however 
that this last two-parameter model is not compatible with the factorisation formula ©, valid at small values of x, 
where only the wavefunctions depend explicitly on Q 2 . 

Moreover, in order to understand to what extend our analysis depends upon the model used for the vector-meson 
wavefunctions, we shall test the dependence of our results with respect to two choices which have already proven to give 
good results. As mentioned before, we will consider the boosted Gaussian (BG) and the light-cone Gaussian (LCG) 
wavefunctions. For completion of the phenomenological discussion, we will also introduce some other combinations of 
wave- functions called BLL and BLB, as explained in the next subsection. 
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TABLE I: % 2 results of the fits to the vector-meson production data. The different wavefunctions indicated (BG, LCG, BLL, 
BLB) are discussed in the text. For the best wavefunction (BLB), the x 2 values are compared with alternative models, that 
feature a t— independent saturation scale. 
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4.077 ± 0.310 
3.754 ± 0.095 


4.472 ± 0.325 
3.724 ± 0.093 


4.258 ± 0.332 
3.708 ± 0.097 


4.041 ± 0.311 
3.713 ± 0.096 
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TABLE II: Parameters obtained from the fits. The parameters of our model do not vary much when using different wavefunc- 
tions. The LCG wavefunction describes better the <j> meson and the BG wavefunction describes better the p meson. This is 
why the BLL or BLB choices give better global descriptions. 

C. Results 

We show in Table[T]the x 2 corresponding to the fits of the different models presented above to the specified data. The 
parameters associated with those fits are given in Table HU In both tables, successive columns show the results for our 
model with t— dependent saturation scale convoluted with BG, LCG, BLL and BLB wavefunctions (we shall discuss 
the wavefunctions labelled BLL and BLB in a moment), followed by the two alternative models without momentum- 
transfer-dependence in the saturation scale (using BLB wavefunction). Those results show that our approach gives 
good results for vector-meson production and this is confirmed by Figures [H and [3] on which we have plotted the 
curves for our fit together with the experimental data points. 

It can be realized from FigsH]and[3]that we obtain a satisfactory description of all data, considering the very small 
number of free parameters. It is of particular relevance for the inclusion of charm in the analysis to obverve that the 
different behaviour of the J cross-section at fixed W from lighter vector mesons is reproduced within the same 
model. Hence the predicted mass dependence (via the rapidity dependence) together with the transfer dependence 
of the saturation scale seem to be both relevant in the successful description of data. Indeed, in our parametrisation, 
the mass dependence is not present in the non perturbative form factor. 

Some technical comments are in order: 

• First we notice that the J/^> differential cross-section always features a bad x 2 , and this is the case for all the 
wavefunctions and models studied. This is probably due to a normalization discrepancy between the HI and 
ZEUS data, which are both included in the fit (in all the cases under consideration, we obtain a good description 
of the ZEUS data and a rather poor description of the HI data) . Fortunately, this concerns only a few points, 
and it does not alter the global x 2 too much. 

• We obtain decent fits with both the BG and LCG wavefunctions though, if we have a closer look at the partial 
values of x 2 , we realise that there is a significant dependence on the choice for the wavefunction. The LCG 
description is significantly better for the 4> meson while the BG description is better in the p meson case. 
Globally, the LCG x 2 1S slightly better than the BG one. 

This inspired the following studies: using simultameously the LCG wavefunction for the <j) meson and the BG 
wavefunction for the p meson. When the LCG wavefunction is used for the J /"J the analysis is labeled BLL, 
and when the BG one is used, it is called BLB. 




FIG. 2: Fit results for the p (HI 0), 4> (ZEUS [H) and J/* (ZEUS Q) elastic cross-sections. 



One sees that these combinations really improve the \ 2 values (this is the case for all the three models under 
study). In addition, the parameters c and B do not depend much on the choice of the wavcfunction. Hence, 
this represents a good global description of the data and shows that the confrontation of our model with the 
experiments is successful. 

Comparing fits performed with our parametrisation (|12p and the corresponding one with a t — independent 
saturation scale and constant slope, one sees that the introduction of Q s (t) significantly improves the fit. This 
significative result, which is consistent with theoretical expectations, opens the way for further experimental 
test of the saturation regime of QCD with exclusive processes. 

Compared to our prametrisation, the parametrisation with a t— independent saturation scale and with B{Q 2 ) = 
B + B' /(Q 2 + My) shows \ 2 values lower in the 0-meson case but higher in the p- and J/'J'-meson case. The 
global x 2 is slightly higher than that of our QCD-inspired model, in which this effective behaviour in Q 2 is 
accounted for by t dependence of the saturation scale. 

In all cases, we have used the BLB choice for the wavefunctions as it gives the best results. By doing this 
we hope to minimise the impact of the wavefunctions in our discussion and focus on the description of the t 
dependence of the scattering amplitude. 

For better comparison with the other parametrisations, we have displayed in Figure U] the curves corresponding 




FIG. 3: Fit results for the p (HI gtj), <j) (ZEUS [p]) and J/* (ZEUS [H, HI [13]) differential cross-section. 



to the three parametrisations of Table [J (our i-dependent parametrisation for the BLB wavefunctions together 
with the two i-independent ones) for p-, (f>- and J/^-meson differential cross-sections for the lower and higher 
available values of Q 2 . 

• The plot on figure [3jc) deserves a special comment. The data are for the differential cross-section at t = 0, 
therefore the corresponding theoretical curves only depend on the parametrisation of the forward scattering 
amplitude, and on the choice for the 0-meson wavefuntion. We have observed a x 2 P er point of 3.157 (resp. 
1.625) for the BG (resp. LCG) wavefunction, suggesting again that the LCG wavefunction might be more 
appropriate to study the <p meson production. 

• From the value of the parameter c = 4 GeV -2 obtained in our approach, we can say that the saturation scale 
starts to increase with momentum transfer when t > c _1 = 0.25 GeV 2 . It is interesting that this scale lies within 
the range of accessible data meaning that we can get insight on the dynamics of saturation by looking at the 
present HERA data. 
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FIG. 4: Comparison of the three parametrisations for differential cross-sections. The left plot shows the p-meson production at 
Q 2 =3.1 and 21.2 GeV 2 , the middle one displays the ci meson at Q 2 =2.4 and 19.7 GeV 2 , and the rightmost corresponds to the 
J/ty meson at Q 2 =3.1 and 16.0 GeV 2 . Continuous lines: t-dependent saturation; fat-dashed lines: t-independent saturation, 
fixed slope; thin-dashed lines: Q 2 -dependent slope. 



W=82 GeV 





10 1 10 2 

Q 2 (GeV 2 ) 

(a) 7*p — ► 7p at fixed W 

FIG. 5: Predictions for the DVCS measurements. The data are from HI [2 
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D. Predictions for DVCS 



Using the overlap function {7} in formula ((6]), we can predict the DVCS cross-section. Those predictions obtained 
with our model and the parameters determined above are presented in Figure They are compared with a few 
available data [25| and the agreement is good. Because of the rather large error bars, these data do not alter the fit 
described previously, when included. Note that contrarily to diffractive vector-meson production, the dipole-model 
description of the DVCS data involves a wavefunction that can be computed from perturbative QED and does not 
require a model. The DVCS predictions of the different models presented in the previous sections differ quite a lot, 
at least in shape, as seen in Figure We hope that forthcoming DVCS data will help to further test our predictions. 



V. CONCLUSIONS AND PERSPECTIVES 



Following the theoretical observation that the geometric scaling property of the BK saturation equation is valid not 
only for inclusive but also for diffractive exclusive processes of DIS, we have proposed an extension of the parametri- 
sation [l8[ (see formula (|12[l ) to non- forward vector- meson production. The theoretical analysis predicts a saturation 
scale proportional to t and a factorisation of the nonperturbative proton form factor. Hence, the geometric scaling 
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FIG. 6: Predictions of the three models studied in this paper for the DVCS diffractive cross-section. For more clarity, we have 
multiplied the differential cross-section by exp(5|i|). The thin lines are related to Q 2 = 4 GeV 2 and W = 71 GeV, while the 
thick ones are obtained with Q 2 = 8 GeV 2 and W = 82 GeV. 



property, i.e. the Q/Q S (Y) scaling, is preserved at non-zero momentum transfer with Q s (t, Y) = Q s (Y)f(t), and f(t) 
interpolating between the soft scale at small t and \fl at higher t. 

By introducing two parameters to feature the predicted behaviour of the saturation scale with t (see (jTUJ) and 
the factorised non perturbative proton form factor (see (fT^f ). we satisfactory describe the available HERA data for 
exclusive diffractive p-, <f>- and J/'F-meson production. This good agreement with the data shows the consistency 
between QCD saturation predictions and measurements in the HERA energy range. We shall be able to further test 
our parametrisation when the preliminary p, <j>, and DVCS data become final. Note also that the cxtcntion (fT2|) to 
non- forward amplitudes could be considered for t = saturation parametrizations other than pip . 

Other successful models of exclusive processes with saturation effects have already been achieved [2fj| . [H, [29| 
confirming the interest of saturation physics for vector-meson production. However we would like to emphasize the 
specific aspect of the present work, which may open a new way to describe exclusive measurements in DIS. Indeed, 
within our approach, we use the momentum transfer q instead of the impact parameter b to parametrise the saturation 
scale. This q-dependence is expected from p-QCD while the nonperturbative dependence in q is factorised, which is 
not the case in impact-parameter space. This is also practically convenient since the data are directly measured as a 
function of t = — q 2 . 
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APPENDIX A: PHOTON AND VECTOR-MESON WAVEFUNCTIONS 



For the self-consistency of the paper, let us give the expressions for the overlap functions we have used in © for 
the different processes. According to equation ([T]), they are the product of a wavefunction for the splitting 7* — * qq 
and a factor accounting for the production of the final state qq — ► 7*, 7, V, appropriately summed over the helicity 
and flavor indices. 

When the final state is a photon (real or virtual), all vertices can be computed from perturbative QED at lowest 
order in the electromagnetic coupling. The results are 



<Mz,r,Q 2 ) = £\ 



*V {z,v,Q 2 ) 



$^(z,r,Q 2 ) 



2 a e N c 
' f 2tt 2 

/ 

E 2 a e N c 

f 

E e /^T {[** + (! - ^QfK^rQ^mfK^rmf) + m 2 f K (rQ f )K (rm f )} , 



4Q 2 z 2 (l-z) 2 K 2 (rQ f ) 



{[z 2 + (1 - zf]Q)Kl(rQ f ) + mjKUrQf)} , 



where e / and m f denote the charge and mass of the quark with flavor / and with Q 2 = z(l — z)Q 2 + m 2 ?. In practice, 
we sum over five flavors. 

When the final state is a vector meson, the light-cone wavefunctions are usually parametrised in terms of an 
additional unknown vertex function for which there exists different models. The expressions to be used in ^ are 
then (for a vector-meson V) 



< y (z,r,Q 2 ) = e fx /^N c 2QK (rQ f ) 



lOte 



M v z(l- z)4> L {r,z) + 6- 



7 



<£l' v (z 7 r 7 Q 2 ) = &f\I^N t 



OeNc 
27T 2 



M V 

2l 



-4>L(r, z) 



{m 2 f K {rQ } )Mr,z) - [z 2 + (1 - z) 2 ]Q fK^rQ f )d r T {r, z)} 



where the constant if stands for an effective charge. It is given in table linl along with the quark and meson masses 
used. Those expressions are very similar to the ones use for photons except for the function (j>L,T that is related to 
the vertex function and depends on the model. The first choice we have investigated is the boosted Gaussian (BG) 
wavefunctions which is a simplified version of the Nemchik, Nikolaev, Predazzi and Zakharov model [l5j with (5 = 1 
and 



i>L,T = N LtT exp 



m)R 2 
iz{l~ z) 



mjR 2 



2z(l- z)r 2 



R 2 



The parameters R and Nl^t are constrained by unitarity of the wavefuntion as well as by the electronic decay widths. 
They are given in table Hill 

An alternative solution is the light-cone Gauss (LCG) wavefunction [lTj]. For those, we take (5 = and 

<j> L = N L exp [-r 2 /(2R 2 L )] , 

cPt = N T z(\ - z) exp [-r 2 /(2i? 2 )] , 

with the parameters also given in table Hill 



Vector-meson 


common parameters 


BG parameter 


s 


LCG parameters 


M v (GeV) 


m f (GeV) 


e/ 


R 2 (GeV" 2 ) 


N L 


iVr 


R\ (GeV 


- 2 ) 


R\ (GeV" 2 ) 


N L 


N T 


P 


0.776 


0.14 


1/V2 


12.9 


0.853 


0.911 


10.4 




21.0 


1.79 


4.47 




1.019 


0.14 


1/3 


11.2 


0.825 


0.919 


9.7 




16.0 


1.41 


4.75 


J/* 


3.097 


1.4 


2/3 


2.3 


0.575 


0.578 


3.0 




6.5 


0.83 


1.23 



TABLE III: Parameters for the vector-meson light-cone wavefunctions. 
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